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We model disorder in graphene by random impurities treated in a coherent-potential approxima- 
tion. Using the analytically solvable Lloyd model for the disorder distribution, we show that the 
temperature dependence of the minimum conductivity as well as the temperature dependence of 
the resistivity at high densities and the density dependence of the respective slopes are consistently 
explained by a temperature dependent disorder strength _T consisting of a constant plus a T-linear 
contribution. This finding suggests that at least two contributions to scattering in graphene are 
important for its transport properties, and that one of the contributions is due to scattering of 
electrons from thermally induced excitations. 

PACS numbers: 72.80.Vp, 72.10.-d 
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I. INTRODUCTION 

Two of the hallmark features of graphene are its linear 
Dirac-cone quasiparticle dispersion and the finite min- 
imum conductivity^ In non-suspended graphene (NSG) 
the minimum conductivity <r m i n , he., the minimal value of 
the dc conductivity er^c with respect to variations of the 
electron density, is weakly dependent on temperature T 
but varies considerably from sample to samplei^ In the 
attempts to approach the ballistic limit of Dirac fermions 
without scattering, realizations of suspended graphene 
(SG) sheets have been prepared, which enable unprece- 
dented electron mobilities^ and show a stronger increase 
of er m i n upon increasing T£& For low electron densities 
when the Fermi energy is slightly off-center of the Dirac 
cone, <7dc decreases with T below a crossover temper- 
ature, while the resistivity increases linearly with T at 
high densities^ 

Free electrons on the half-filled honeycomb lattice con- 
stitute a perfect conductor with a non-zero Drude weight 
at finite temperature and a finite dc conductivity <jq — 
ire 2 /2h at zero temperature^— This is a direct conse- 
quence of the Dirac-cone structure of the electronic dis- 
persion. For finite temperatures also the optical conduc- 
tivity is of order Co in the visible frequency range. This 
theoretical result, neglecting disorder effects, is indeed 
observed in optical absorption experiments on charge- 
neutral graphene ^ 

Graphene samples are not pristine, however ^ 2 : 13 
Scanning-tunneling microscopy images of NSG samples 
show inhomogeneous patterns^ their origin was traced 
to the presence of impurities at the substrate-graphene 
interfaced Impurities may nucleate electron- or hole- 
rich puddles^ which obscure the intrinsic Dirac fermion 
physics of pristine graphene. In suspended graphene 
(SG) scattering can occur due to microscopic corruga- 
tions of the otherwise unstable two-dimensional crys- 
tal, so-called ripples, which were observed by transmis- 
sion electron microscopyii and theoretically analyzed as 



one possible source for electron scatterin g 18 : 19 or charge 
inhomogeneity^ Impurity effects are quantitatively less 
important for the optical conductivity, but play a ma- 
jor role in determining the dc transport properties of 
graphene; the latter are the topic of this paper. 

o"dc of graphene is minimal at the charge neutrality 
point, which corresponds to a honeycomb lattice at half- 
filling, i.e., with one electron per lattice site. Upon 
applying a gate voltage of either sign, Ode increases^ 
hence the name "minimum conductivity" . Disorder- 
induced charge-density modulations imply a spatially 
varying chemical potential and thereby conceal clean 
Dirac fermion physics. Also the dc transport measure- 
ment process itself may introduce a bias, e.g., due to a 
charge transfer at metal contacts^ 

Theoretical work on transport in graphene com- 
prises studies of charged-impurity scattering, as reviewed 
in Ref. |22, different sources of disorder ,2 s - and also 
the crossover between low- and high-density regimes* 2 ^ 
More recent studies focused on ballistic^ or diffusive 
transport*^ and also the effects of finite-range scatter- 
ing at finite densities ^ Unresolved problems remain in 
particular in the low-density regime, which is relevant for 
the minimum conductivity at zero bias. 

Several predictions exist for a minimum conductivity 
of 4e 2 /nfi in the absence of disorder i22r— This limiting 
value at zero temperature is obtained if the dc limit is 
taken first and the zero-disorder limit afterwards i 9 . Ex- 
periments on both NSG and SG samplesi^ with non- 
universal values of the minimum conductivity were re- 
ported, with the trend that the T dependence of the 
minimum conductivity is enhanced in clean SG samples 
as opposed to dirty SG£ or NSG samples^ In fact, the 
minimum conductivity increases with increasing temper- 
ature, i.e., as in a semiconductor. In contrast, at suffi- 
ciently large gate voltages a metallic T dependence of the 
conductivity is observed, with a resistivity increasing lin- 
early with T and a slope that decreases upon an increase 
in the gate voltage^ 
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FIG. 1. Transport at the charge neutrality point in pristine 
graphene. (a) At T = the valence band (states below the 
chemical potential /i), shown here for a single Dirac cone, is 
filled (indicated by the dark blue shading) and the conduc- 
tion band is empty. At T > electrons are thermally excited 
to the conduction band (lighter blue shading), (b) Contri- 
butions to the conductivity: Only interband transitions are 
allowed at T = due to Fermi blocking. Intraband transitions 
contribute to dc transport at finite temperatures, (c) Dynam- 
ical conductivity of pristine graphene within the Dirac-cone 
approximation at T = and T > 0. At T = interband tran- 
sitions lead to a universal finite conductivity of crn. = 7re 2 /2/i. 
For T > a Drude peak emerges due to the intraband transi- 
tions with a Drude weight D oc T. For visual frequencies OMs 
the optical conductivity still is of order crn. 



Here we evaluate the Kubo formula for the dc conduc- 
tivity of electrons with a linear Dirac cone dispersion. 
Disorder effects are included by a random chemical po- 
tential, which is treated within the coherent-potential ap- 
proximation (CPA). 33 ' 34 The associated disorder energy 
scale r may itself depend on temperature. We specifi- 
cally investigate the case of a Lorentzian disorder distri- 
bution of width r ("Lloyd model"), for which the Kubo 
formula can be evaluated exactly within CPA. As a re- 
sult <7dc at half-filling depends only on the type of dis- 
order distribution and the dimensionless ratio T '/ T '. For 
T / r S> 1 the minimum conductivity increases linearly 
with T/r. For a temperature dependent r = i~b + a{T 
the minimum conductivity thus saturates at high tem- 
peratures. With this simple ansatz and a choice of typ- 
ical meV energy scales for r, the T dependence of <7d c 
changes from semiconducting at half-filling to metallic at 
sufficiently large band filling. At intermediate densities 
Cdc evolves from metallic to semiconducting behavior in 
the temperature range between and 200 K. Moreover, 
the experimentally observed linearly increasing resistiv- 
ity at high temperatures in the metallic regime as well 
as the decreasing slope upon increasing the density are 
reproduced in this ansatz. 



II. MODEL AND METHOD 

An infinite sheet of pristine graphene is modeled by a 
tight-binding Hamiltonian with an effective next nearest 
neighbor hopping on a honeycomb lattice without impu- 
rity scattering and electron-electron interactions. In the 
absence of current- vertex corrections the dc conductivity 
follows from 
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(1) 
where f x = 1/(1 + exp(a;/T)) is the Fermi-Dirac 
distribution function (with ks = 1) and /5(e) = 
L^ 1 J2k(® 6k I 'dk x ) 2 S(e— efe); L is the number of unit cells 
of the lattice, /i is the chemical potential which vanishes 
at half-filling. In Eq. (pj a prcfactor of 4 has been incor- 
porated; it accounts for the spin and valley degeneracies 
of graphene. 

For free electrons the spectral functions simply reduce 
to A c {y) = 5{e—v). We use the Dirac cone approximation 



p(e) = h\e\/2w for 



e < e n 



where e max is a cutoff en- 



ergy chosen as the half-bandwidth of graphene. Indeed, 
the Dirac cone approximation for /5(e) gives the correct 
result for <7dc and serves as a good approximation even in 
the visual frequency ranged where the band dispersion 
leads to only weak quadratic corrections to (j(lo) at low 
frequencies. 

The term in Eq. ([!]) which involves A c (v) 2 in the in- 
tegrand leads to the usual intraband conductivity as in 
single-band models. It gives rise to a Drude-like con- 
tribution, hence an infinite dc conductivity in a perfect 
conductor. Also in the presence of electron-electron in- 
teractions this expression for the intraband conductivity 
remains correct, if the self-energy S{v) is local and A e (is) 
= — Im(e — v — Z , (^))~ 1 /7Ti 35 i 36 The second term in Eq. 
([T]), involving A e {y)A_ e (y), describes excitations with a 
particle at energy e and a hole at — e and accounts for 
interband transitions. This contribution accounts for the 
visual transparency of graphene in the dc limit of the 
optical conductivity, <r = ire 2 /2h (see Fig. Hfc)Ai 

The presence of disorder complicates the situation con- 
siderably. Discrete translational invariance is broken, 
rendering microscopic theoretical approaches much more 
difficult than in the homogeneous case. One standard 
approach is the Anderson model 3 -^ with local potential 
impurities, 



h = h q + Y, Km 



(2) 



where Hq is the tight-binding Hamiltonian for the clean 
system, rii is the local density operator on site i of the 
lattice and Vi is a random variable determined from a 
probability distribution P(Vi). Here we do not aim at a 
full microscopic description of disorder, e.g., in the spirit 
of a self-consistent diagrammatic treatment of impurity 
scattering effects^ and recall that weak localization is 



suppressed by long-range scattering in graphene^2r— 
Instead we apply the coherent-potential approximation 
(CPA) 33 ' 34 ' 42 to determine an effective random medium 
described by a local self-energy S(u>), which is deter- 
mined by a self-consistent solution of the CPA equations 
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In Eq. (O G Q (z) = /dw Pdos(u)/(z - u) is the lo- 
cal Green function of the clean system described by Hq, 
Q{oj) is a dynamical Weiss field and D[z] is the Hilbcrt 
transform with respect to the disorder distribution func- 
tion P(Vi). The CPA expression for the conductivity^ 
agrees with the Kubo formula Eq. (JTJ with A e (v) = 
-Im(e-i/- S{v))- 1 /tt. 

The CPA equations ([3]) can be solved, at least numer- 
ically, for an arbitrary disorder distribution. In order to 
keep the subsequent analysis as simple and transparent 
as possible, we focus on the specific case of a Lorentzian 
disorder distribution of width I\ 
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This is the so-called Lloyd model, for which the CPA 
equations are exactly solvable using D[z] = (z + i-T) -1 , 
which yields S(v) = —IT. Hence we obtain for the Lloyd 
model 
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as the input quantity for Eq. (JTJ. A e (y) is of the form 
a{{e — v)/r)/r i implying that the dc conductivity is only 
a function of the ratio T/T for fx = and of [xjF for 
T = 0; the latter holds only, if \i <g e max , which is fulfilled 
in the experiments cited above. 

The scaling behavior of <7 m i n has two reasons: (a) the 
cutoff energy e max can be replaced by infinity in the e- 
intcgral in Eq. (TTJ and thus docs not appear as an addi- 
tional energy scale, and (b) for dimensional reasons the 
dc conductivity is universal in the sense that it does not 
depend on the hopping matrix element t of the underly- 
ing two-dimensional tight-binding Hamiltonian and, con- 
sequently, not on the Fermi velocity vf- Both reasons 
are directly related to the linearity of the dispersion in 
graphene up to energies much larger than the relevant 
temperatures. For finite densities this universality no 
longer holds, since the density variations are determined 
by the chemical potential which thereby depends on the 
hopping matrix element t. 

Experimentally it is the gate voltage which controls the 
electronic density n, measured relative to half-filling. For 
given temperature T, disorder strength r, and chemical 
potential \i the density is given by 



dw PDOSM (fu-fi - fu) , 



(6) 



PDOsM 



V3nt 2 A u ./_ 



dc \e\A e (uj) (7) 



is the density of states (summed over both spin projec- 
tions) for the disordered system, t = 2.7 cV the hop- 
ping matrix element of the tight-binding model, A u = 
3V3a 2 /2 the size of the unit cell, and a = 1.42 x 10~ 10 
m the interatomic distance on the honeycomb latticed 



III. RESULTS 

We first keep the disorder strength _T fixed and discuss 
basic properties of the minimum conductivity and the 
conductivity at finite chemical potential and zero tem- 
perature. In a second step we evaluate the density depen- 
dence of the conductivity for a typical disorder strength 
(on the order of meV— ^) and for temperatures in the 
range from K to 200 K. Especially we consider the 
temperature dependence of the resistivity crj c at fixed 
densities and show that a temperature-independent r at 
high densities is insufficient to explain the experimentally 
observed linear T dependence of the resistivity. Adding a 
phenomenological linear T dependent contribution to r, 
the experimental observation is matched by our ansatz. 
Moreover, the observed density dependence of the slope 
in the T-linear regime of the resistivity follows naturally 
without further assumptions. For a selected T depen- 
dent r = Tq + a,\T , the minimum conductivity increases 
with temperature but with a decreasing slope; er m i n sat- 
urates at high temperatures closely similar to the exper- 
iments. Finally we show the T dependence of the con- 
ductivity with a density dependent crossover from metal- 
lic (dtTdc/dT < 0) at low T to semiconducting behavior 
(dcr dc /dT > 0) at high T. 

The temperature dependence of the minimum conduc- 
tivity <7 m i n for the Lloyd model with disorder strength 
r is shown in the top panel of Fig. [2] For T/T -> 0, 
Cmin tends to the limiting value 4e 2 /irh, which coincides 
with the clean limit discussed in Refs. I28H33 . However, 
Ae 2 /nh should not be considered a universal value, but 
rather a particular result of the Lloyd model. Other 
non-Lorentzian disorder distributions arc likely to lead 
to other values of the minimum conductivity. At high 
temperatures cr m in increases linearly with T jT for T jT 
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To understand the physical processes involved we dis- 
cuss the relevant contributions to <7 m i n - First we note 
that the clean case at zero temperature is not recovered 
by our theory for <j^ c . However, this is not a shortcom- 
ing but rather a generic feature of the conductivity as 
a function of disorder strength, temperature, and fre- 
quency. We recall that at zero temperature intraband 
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FIG. 2. DC conductivity in units of e 2 /h for the Lloyd 
model in CPA. Top panel: Minimum conductivity (/i = 0) as 
a function of T/T. Bottom panel: DC conductivity at T = 
as a function of p/T. 



excitations are prohibited. Thus only the interband ex- 
citations are responsible for Re u(a; ^ 0) = <jq^ (see Fig. 
[T]), i.e. when the dc limit is taken after the limits of zero 
temperature and zero disorder strength. The theory pre- 
sented here instead aims at describing dc measurements, 
for which the dc limit must be taken first. In the latter 
case, both interband and intraband excitations are rele- 
vant and both contribute equally (2e 2 /wh for the Lloyd 
model) to the T/T —t limit. The discrepancy between 
<7o and c m ; n (T — > 0) may also be understood by noting 
that for do the largest energy scale in the system is the 
frequency (taken to zero last), while the largest energy 
scale for er min is the disorder strength 7". 

At finite temperatures thermally excited particles in 
the conduction band render intraband particle-hole exci- 
tations possible, leading to a non-zero Drude weight and 
thus an infinite <7dc (for r — > and therefore T '/ T — > 
oo), while interband low-energy excitations are blocked 
by thermally occupied states in the conduction band. 



Viewed as a function of temperature at fixed r the con- 
ductivity is semiconducting, i.e., d(7dc/dT > 0. 

At finite chemical potentials the interband excitations 
become less important, and the behavior is determined 
mostly by intraband excitations. The situation for higher 
densities thus resembles more and more the case of a 
single partially filled band, where metallic behavior sets 
in for sufficiently low temperatures, i.e., Aa dc /dT < 0. 

The dc conductivity at T = as a function of p/T is 
shown in the lower panel of Fig. [2] It tends to the limiting 
value 4e 2 /7r/i for \p/F\ — > and increases linearly for 
\p/T\ 3> 1 but well below the cutoff e m ax/^- This linear 
dependence on the chemical potential for large \p/r\ is 
analytically obtained from the Kubo formula, taking into 
account intraband excitations only, 
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where the last asymptotic expression is valid for \p/r\ 3> 
1. 

For fixed value of the disorder strength r = 1 meV we 
show in Fig. [3] the dc conductivity and its inverse, the 
resistivity p, as a function of density for selected tem- 
peratures. Since the density depends quadratically on 
the chemical potential for p/T ^> I and Ode depends 
linearly on |/i/-T| in this limit at zero temperature, the 
low-temperature conductivity increases like y/n at high 
densities. Experimentally a sublinear density dependence 
of the conductivity was also reported in Ref. y. 

In Ref. [fj a linear increase of the resistivity as a func- 
tion of temperature was observed for high densities at 
elevated temperatures. For a temperature independent 
r, <7dc = \/J>/r\ e /h for large \/j,/r\, and the tempera- 
ture dependence of the chemical potential at fixed den- 
sities follows fi = fj,(T = 0) + 0{T 2 ). Hence, a dc (T) 
= a dc (T = 0) + 0(T 2 ) for fixed densities and thus p = 
p(T = 0) + 0(T 2 ). A linear temperature dependence of p 
therefore requires a temperature dependent _T within our 
ansatz. In fact, a linearly increasing resistivity at high 
densities naturally follows from _T = Jq + a\T , 



r h /q h ct\T h 



p e* 



p e' 
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(10) 



p thereby depends not only explicitly on temperature, 
but also implicitly via the T dependent r. This implicit 
T dependence is, however, negligible for large fillings, 
when also the explicit T dependence is very weak since 
it scales like temperature over Fermi energy 

In the following we adopt the T dependent disorder 
strength r(T) = _Tb + a± T and fix the parameters To = 
1 meV and a\ = 3.25 meV/200 K such that the tempera- 
ture dependence of the minimum conductivity (see lower 
panel of Fig.JSJ) approximately matches the experimental 
data of Ref. [fl The density dependence of the resistivity 
for the selected T dependent disorder strength is shown 
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FIG. 3. Upper panel: Conductivity as a function of density 
at different temperatures for a temperature independent F. 
Lower panel: The data for the resistivity p — a^ . 
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FIG. 4. Upper panel: Resistivity as a function of density for 
different temperatures. Lower panel: Resistivity increase Ap 
= p{T) — p(10 K) above 10 K as a function of temperature 
for different densities. The dashed lines are linear fits to the 
data points between 100 K and 200 K. 



in the upper panel of Fig. 01 For moderate temperatures 
below 200 K the temperature dependence of the resis- 
tivity is presented in the lower panel of Fig. 2) Eq. (fTU)) 
implies that the slope in the linear regime is proportional 
to ai/fx, and since \i ex y/n the slope decreases like 1/y/n. 
For high temperatures well above 200 K and sufficiently 
large densities, or below 200 K for moderate densities, 
there is a deviation from linear behavior, and the resis- 
tivity decreases again due to interband excitations. 

The crossover from metallic behavior at finite densities 
and low temperatures to semiconducting behavior at el- 
evated temperatures is shown in the upper panel of Fig. 
[5J The crossover temperature vanishes at zero density (/i 
= 0), since da m i n /dT > for all temperatures, and also 
increases with increasing density. In fact, for the selected 
temperature dependence of r the conductivity is metal- 
lic below 200 K for densities larger than 8 x 10 10 /cm 2 . 
The temperature dependence of the minimum conduc- 
tivity for the same T dependent r is shown in the lower 



panel of Fig. [5] Here a sublinear T dependence of cr m i n 
is observed for elevated temperatures. Indeed, the curva- 
ture of Cdc changes sign at an intermediate temperature 
depending on the relative sizes of To and a±. An increas- 
ing (7 m i n as a function of temperature with a sublinear 
behavior at elevated T, yet below 200 K, is similarly ob- 
served in experiments^^ 



IV. SUMMARY AND DISCUSSION 

We have presented a phenomenological theory for 
the temperature dependence of the dc conductivity of 
graphene at zero and finite particle densities including 
potential disorder in a coherent-potential approximation 
(CPA). Specifically we have chosen a Lorentzian disorder 
distribution ("Lloyd model"), for which the CPA equa- 
tions are exactly solvable. This approach recovers well- 
established limits in the clean case and at the same time 
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FIG. 5. Upper panel: Conductivity as a function of tem- 
perature for different densities with a T dependent F in a 
double-logarithmic scale. The density is measured with re- 
spect to half-filling. Lower panel: Temperature dependence 
of the minimum conductivity. 



provides a phenomcnological context for the remarkable 
transport properties of graphene in the presence of impu- 



rity scattering. For the Lloyd model the minimum con- 
ductivity is 4e 2 /irh, which coincides with previous pre- 
dictions for the dc limit in the clean system provided that 
the zero frequency limit is taken before the clean limit at 
zero temperature. At finite temperatures, the enhanced 
T dependence of the minimum conductivity in cleaner SG 
samples is explained, and we find er m m ex T/T for T/T 
^$> 1. As a consequence we expect a very steep increase 
of "min with temperature in even cleaner samples. More- 
over we have shown that the T linear resistivity at high 
densities and the density dependence of its slope follow 
naturally from a temperature dependent r = Jo + a\T. 
This phenomenologically determined T dependence of r 
suggests the existence of at least two sources for scatter- 
ing in suspended graphene. The constant Jo points to 
static potential disorder, whereas the T-lincar part u\T 
may arise from scattering off a thermally excited per- 
turbation. One obvious possibility are thermally excited 
ripples, since even the linear T dependence of the scatter- 
ing rate could be explained within the ripple scenario^. 

Here we have investigated the role of disorder, as de- 
scribed by an Anderson impurity model with a phe- 
nomcnological disorder strength, as a source for scat- 
tering in graphene. As pointed out in Ref. [2l], it is 
important to understand which additional extrinsic ef- 
fects may mask the intrinsic properties of graphene, es- 
pecially the sensitive Dirac fermion physics at the neu- 
trality point. Possible extrinsic perturbations are con- 
tact resistances, spurious chemical doping into the con- 
tact regions, or macroscopic charge inhomogeneity on 
length scales comparable to the sample size. Such ef- 
fects need to be incorporated in order to understand the 
unusual transport properties of graphene in particular at 
the charge-neutrality point. Also improved doping tech- 
niques using organic molecules^ may help to unveil the 
intrinsic transport properties of grapheme. Further theo- 
retical and experimental activity should clarify these as- 
pects and the promising prospects of graphene as a basis 
of future electronic devices. 
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